Thermodynamics of the superfiuid dilute Bose gas with disorder. 
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We generalize the Beliaev- Popov diagrammatic technique for the problem of interacting dilute Bose 
gas with weak disorder. Averaging over disorder is implemented by the replica method. Low energy 
asymptotic form of the Green function confirms that the low energy excitations of the superfiuid dirty 
Boson system are sound waves with velocity renormalized by the disorder and additional dissipation 
due to the impurity scattering. We find the thermodynamic potential and the superfiuid density 
at any temperature below the superfiuid transition temperature and derive the phase diagram in 
temperature vs. disorder plane. 



Superfluidity in random environments enjoys a long 
standing yet intense attention. The effect of disorder 
on the behavior of systems possessing long-range cor- 
relations is central to contemporary condensed matter 
physics, and superfiuid Bose gas offers an exemplarily 
unique and accessible tool for both experimental and the- 
oretical researches. One of the fascinating properties of 
such systems is their ability to maintain superfluidity (i.e. 
long range correlations) even in the strongly disordered 
environment. He 4 , for example, remains superfiuid when 
absorbed in porous media [Q. The problem of influence 
of disorder on superfluidity (and on its close analog - 
superconductivity) has been under extensive theoretical 
attack (see seminal works ) and remarkable progress 
in qualitative understanding of disordered Bose systems 
was achieved. Recent papers [|],|| discussed a continuum 
model of the dilute interacting Bose gas in a random po- 
tential. The advantages of this model are that (i) it is 
microscopically related to the original problem and (ii) it 
is very well understood in the clean limit. The proposed 
model describes, in particular, the quasiparticle dissipa- 
tion and depletion of superfluidity at zero temperature 
and marked an important step towards quantitative de- 
scription of disordered Bose systems. 

In this Letter, building on the model of Refs. 
we develop a systematic diagrammatic perturbation the- 
ory for the dilute Bose gas with weak disorder at finite 
temperatures below the superfiuid transition tempera- 
ture T s . We obtain disorder corrections to the thermo- 
dynamic potential which completely determine thermo- 
dynamic properties of the superfiuid system. We derive 
for the first time the disorder-induced shift of T s result- 
ing from disorder scattering of quasiparticles with energy 
e ~ T. We find that the superfiuid density decreases 
monotonically with the temperature. This completely 
agrees with the experimental data, while being in some 
contradiction with the theoretical result of Ref. Q where 
a non-monotonic temperature dependence of superfiuid 
density dependence was claimed. In the limit T — > our 
theory reproduces all the results of Refs. ||||. 

The model. The starting point of our model is the 
Lagrangian density: 



C = -tp*(-Vl/2m — fj, + u(r) + d T )ip - g <p*<p*(p(p, (1) 

where (p — ip(r, r) is the field representing Bose parti- 
cles, r is the real space coordinate, r is the Matsubara 
time and u(r) is the disorder potential. As usual we con- 
sider a soft interaction potential g(r) and use the Born 
approximation g = 4ir\/m to relate the interaction con- 
stant g = J g{r)d d r to the scattering length A. || Taking 
Gaussian (S-correlated disorder (u(ri)u(r2)) = v8{r\— T2) 
we derive the effective replicated action in a form: 



S = 



p 



<f a (p)(k 2 /2m -fi- iio)^ a {p) + V + V d) (2) 



with p = (fc, to) and the interaction parts Vj and Vd 



Vi = 



fJ 



2f3V 

Vd = ~2~V ^ l P*a( k l, UJ l) l Pa(k3,LJl)<Pp{k2,UJ2)<Pf3(k4,U2) 

where a, (3 = 1, ...,n are the replica indices, k = v 2 and 
the conservation of total momentum k\ + k% = ka + fej 
(in Vi and V^) and of total 'energy' 0^1+^2 = + 
W4 (in Vi ) is assumed. The corresponding vertices are 
presented in Fig. la. Below T s we separate the condensate 
contribution by shifting the fields ip a — > a^J [5V 5^,0^,0 + 
(f' a and define the Green functions of the fields <p' a by 

gmp) = (va(py;(p)>,-p^(p) = <^(pH(-p)) • ( 3 ) 

Defining also the functions G a p(p) — (tp 1 * (— p)(p'p{— p)), 
F a p = {fa (~ P)f'p (p)) we introduce the matrix Green 
function and the corresponding matrix self energy 



Qa,j3 



G a 



(3 fa/3 



B a /3 A aj3 



that are related by the Dyson equation 

G^ 1 = (p 2 /2m - fi - iLUT 3 )S a/ 3 + E, 



(4) 



(5) 



where T3 is the Pauli matrix. The condensate density 
is uniformly distributed in the replica space, therefore 
Green function Q and self energy E can be presented as 



1 



Gaf3 = Ql S a f3 + Ql Raft, S Q/ 3 — Si 5 a p + S 2 R a (6) 

where R a p is the matrix with all elements equal to 1. 
From the n — ► limit of Eq.(||) we obtain that the replica 
diagonal part of the Green function is determined by the 
replica diagonal part of the self energy 

Gi 1 = (p 2 /2m - (j, - iujT 3 ) + Si, Qi = -£?iS 2 (?i, (7) 

and the poles of the Green function Q are determined 
by the poles of the Green function G\ . Due to the Gold- 
stone theorem, below the condensation temperature the 
function G\ has a pole at p = 0, and from (0) we find 



Ai(0)-Bi(0)=|i. 



(8) 
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FIG. 1. Interaction and disorder vertices (a) and first order 
contributions to self energies (b). Dashed lines in (b) repre- 
sent the condensate. 
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The contribution ri\ represents the quasiparticle density 
due to quantum fluctuations, in the leading order it coin- 
cides with the well known answer for the pure case. The 
contribution 712 represents the density of nonuniform part 
of the condensate. For the theory to be applicable both 
n\ and ni should be much less than the total density 



A' = Xn^ < 1, k' = Km 2 /(8n%V\n) < 1. 



(15) 



The first relation is the usual low-density condition, the 
second one insures that the uniform part of the con- 
densate is not strongly affected by disorder. To re- 
late the condensate density with the chemical potential 
we need to improve the leading order result ga 2 = fi 
considering next order corrections to (0). The sec- 

(2) (2) 

ond order corrections A\ , B\ contain the contributions 
A^ , B^ 2 ''^ from the quasiparticle interactions and the 
disorder contributions Aj 2 ' , B^ d ' linear in k. Correc- 
tions A^ 2 ' l \ B^^ presented in Fig 2a exactly coincide 
with ones studied in for the pure Bose gas. The disor- 
der corrections A^ 2 ' d \ B^ 2 '^ presented in Fig 2b contain 
all the diagrams that are (i) linear in disorder coupling k 
and (ii) have a similar structure with ones shown on Fig 
2a. Here we present the answers for their linear combi- 



nations S 



(2,d) _ ,(2,d) 



= A\*' a> ± B x 



(2,d) 



Our diagrammatic technique parallels that by Beliaev, 
the difference being that we include the disorder vertex 
Vd along with the interaction vertex Vj . The correspond- 
ing diagrams are shown in Fig lb. 



A 1 = 2ga 2 , B 1 =ga 2 , A 2 



Bn 



(9) 



From Eq.(^) we obtain ga 2 = fx and using Eq.(@) for the 
Green functions we have 



fc 2 /2m + /i + iu> 
e 2 (k)+uj 2 ' 

G2(p) = F 2 {p) 



Fi( P ) = 



-fi 



e 2 (fc) W 
nfta 2 



(fc 2 /2m + 2^) 2 



(10) 
(11) 



Re S (2 ' d) (q, w) = gn 2 - — ^ 



fc 4 G'_{k + q,Lo) 



V (fc 2 + 4^m) 2 



(16) 



Re Y/ 2,d \q,uj) — igri2 



k ^ (fc 2 - 8/im) 2 G'(k + q, lu) 



V 



E 



(fc 2 + 4pm) 2 



where G' ± (p) = Re{G x (p) ± Fi(p)). The self energy B is 
real while A is complex with 



nujk 2 /V 8fim — k 2 
e 2 (k + q) +uj 2 (fc 2 + 4^m) 



(17) 



Using Eq.jl^) along with the well known result for 
S^ 2 ''' 1 from Eq.(||) we obtain 



where s(k) = a/ {k 2 /2m) 2 + /ik 2 /m. We see that the 
spectrum of quasiparticles is not affected by disorder in 
the leading order. The Bose gas density is given by 



n = a + ni + r«2 



(12) 



with 



n = 5 E G i (p)> " 2 = 5 E G2 ( fc )- ( 13 ) 



fc,w=0 



Zero temperature. At zero temperature in three dimen- 
sions using Eqs.(|To|, [ll] ) we obtain 



2 , 40 /nA 3 , k v-^ 2mk 2 

» = ga (1 + —J ) +5 n 2 --^ 



V (fc 2 + 4^m) 5 



Combining this equation with Eq.(|l2|) we obtain the re- 
lation between the density and chemical potential 



32 InX 3 . k v-^ 2m k 2 
u = gn[l + —) - — 



3 V ir ' V ^ (fc 2 + 4/xm) 2 ' 



(18) 



Using the thermodynamic relation p = dE/dN along 
with the leading order relation gn = pL we obtain the 
energy 
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FIG. 2. Second order contributions to the self energies 
Ai,Bi due to interaction (a), and due to disorder (b). The 
solid lines represent the Green functions Gi , Fi the crossed 
solid lines represent the Green functions G2 , F2 . The diagrams 
in part (c) represent two contributions to the normal density 
n„ and nl 2 '. 




2nm 



V 



(19) 



which agrees with ||. The integrals over k in Eqs.(|l8| 
, [l9| ) are ultraviolet divergent, it is the consequence of 
the white noise assumption for the disorder correlation 
function. This divergence is not relevant for the low 
energy physics since it could be absorbed in the renor- 
malization of energy and chemical potential: E — > E + 
uk J2k 2m/ k 2 , (i + k J2k 2m/Vk 2 . 

The superfluid density n s can be found from the nor- 
mal density n n which is determined by the transverse 
current-current correlator n s = n — n n . In the leading 
order n n is given by the diagrams shown in Fig. 2c 



,(2) 



(20) 



where n„ is the normal density of the clean system 

n " X) = 3^vJ2p k2 ( G i(p) ~ f i(p)) which after summa- 
tion over Matsubara frequencies may be written as 



1 



12mTV 



E 



sinh 2 [e(/c)/2T] 



(2) 

and n n is the disorder correction 
T 



,(2) - 



\ " 



. ]inV ^2^[Gi(p) - fi(p)]G 3 (p) = -n 2 . 



Thus at zero temperature n„ = and superfluid density 
becomes n s = n — 4rt2/3 in agreement with Refs. Pt PI- 

The second order corrections to the self energies ([16|- 
[Tt| ) can be calculated explicitly for small q, w leading to 
the following low energy retarded Green function 



(21) 



(22) 



Gf(fc,c) = - 



2ickY{k) 1 



lu>0 (23) 



where the sound velocity c is related to the sound velocity 



of the clean system cq by c 



(l+5n 2 /3n), andT( g ) 



k g 4 /24m 2 c 3 7r is the dissipation of quasiparticles due to 
disorder scattering. The dissipation due to quasiparticle 
(phonon) scattering is known to be of a higher power of 
q: T^ pfl \q) ~ q 5 . The results for sound velocity c and 
quasiparticle dissipation T(q) are in agreement with 

Bose condensation temperature. Now we turn to finite 
temperatures. At temperatures above the condensation 
temperature T c the self energy is given by the first dia- 
grams of Figs. 2a, b. 



A, =2 5 n--^G(fc 



(24) 



and Ai,B\ B2 = 0. Taking the sum over k in ( p4| ) we 
obtain the Green function 



Gri(p) = A-- M -fc,. 



(2m)* /=r 



IUJ, 



(25) 



where the = /i — 2gn + ftj^fc 2m/Vk 2 . The density at 
T = T c can be easily obtained from Eq.(p5|) taking into 
account that at this temperature Jl becomes zero 



n = Ca (mT c / 2?r)5 + K T c m 3 /47r 2 . 



(26) 



Solving this equation for T c we find the shift of the con- 
densation temperature due to disorder: 



T c = Tj°) (l - KT^m 3 /67r 2 n 



(27) 



where T c (0) = 27r(n/Ca)t/ m is the Bose condensation 
temperature of the ideal gas with £3 « 2.612. At fi- 
nite 3 only microscopic amount of particles may con- 
dense into a local potential well, this effect leads to the 
smearing of the condensation transition making T c to be 
the crossover temperature between the normal phase and 
a phase where bosons are locally condensed. The true 
phase transition takes place when the chemical potential 
reaches the mobility edge it may be also obtained 
from the condition n s = 0. (see below) 

Thermodynamics at T ~ T s . The self energies at T ~ 
T s are still given by the diagrams presented in Figs. 1 
and 2, but the first diagram of Fig. la should be included 
already in the first order approximation since the density 
of quasiparticle excitations n\ at T ~ T s is of the order 
of total density. This diagram results only in the shift of 
the chemical potential fi — > fx — 2gn\ — /I, and the Green 
functions are still given by the Eqs.([l0 11) but with /1 — > 
ft. The density in the leading approximation is still given 
by Eqs. jT^Jl^ ), but now ni is not a small correction and 
should be therefore calculated with a higher accuracy. 
The main contribution to ni comes from the energies 
e ~ T » ji. At these energies the Green function is given 
by Eq.(]25|) and, thus, the disorder correction to n\ is the 
same as in Eq. (p6|) , i.e. kTw? / 'An 2 . Taking into account 
this correction and using Eq.(p^) we obtain 



3 



Tli 



2tt 



4tt 2 



n 2 



2 2 

nam 2 



(28) 



where T' = T/Tc°\ To relate the chemical potential with 
the condensate density we need to consider the next order 
corrections to the leading order result ga 2 — p, following 
from Eq.(|). Using found in § along with Eq.© 

we obtain 



n 2 



E(n/g)2mk 2 
(k 2 + Ajlm) 2 



3fi 2 m 2 T Km 3 T 
2^ + 2n 2 ' 



and combining this equation with Eq. (E8h we get an equa- 
tion relating density and chemical potential 

n = fi/g - nT'i + p,$m$T/w + An (d) , (29) 

where Au^' is the disorder contribution 

K syy 'lilll,' 2 l> III''' T 

k 



A n {d} = — V 

gV ^ (k 2 + Afim) 2 



An 2 



(30) 



Using the relation N = —dfl/d/j, we eventually obtain the 
disorder correction to the thermodynamic potential: 



v 



K ^ — 



2/im 



4tt 2 ' 



(31) 



Superfluid density at T 



,(2) 



The disorder contribu- 
tion to the normal density tin' at T ~ T s is related to 
ri2 through Eq.(p2|) with n 2 defined by Eq.(^8|) that takes 
into account the chemical potential shift fj, — > p,. Deriving 
the contribution n.„ one needs to consider that accord- 
ing to Eq.(p5|) the spectrum of quasiparticles is affected 
by the disorder at energies e ~ T s , that results in 



(32) 



where tin^ is the normal density of the clean system de- 
fined by Eq.( ^l|) with the spectrum e(k) = {k 2 /2m) 2 + 
p,k 2 /m. Introducing dimensionless condensate density 
a = — T r 2 we write the superfluid density as 

n s /n = n<f l) /n - 4kV/3 - Ak'VttX T'/fJ 3 , (33) 



where n 



(d) 



(d) 
nl = n 



is the superfluid density of the clean system 
rin L \ The dependence of n s on temperature 
for different amounts of disorder is presented in Fig 3. 
Taking n s = in Eq.(^) we obtain the dependence of 
the superfluid transition temperature T s on disorder as 
shown on the phase diagram presented in Fig. 3. The 
dotted line in this diagram represents the condensation 
crossover temperature determined by Eq. (p7|) . 

In conclusion, we have developed a regular diagram- 
matic approach that enables a quantitative description 
of thermodynamics of superfluid dilute Bose gas in ran- 
dom environment at finite temperatures. We have found 
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FIG. 3. The temperature-disorder phase diagram result- 
ing from Eqs.^Q for A' = 0.03 2 . The dashed part of the 
boundary of the superfluid phase corresponding to the region 
K ~ 1 should be understood as an extrapolation. The insert 
shows the superfluid density n s /n dependence on temperature 
T' = T/T c (0) for different amounts of disorder: k' = 0, 0.1, 0.2 
(top to bottom). 



disorder corrections to condensation temperature, ther- 
modynamic potential, and the superfluid density. Our 
results agree favorably with the experimental findings. 
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